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£H , Abstract 

£ : 

By efforts of several reserchers, it is recently established that the dynamical behavior of 
^ ■ the cosmological perturbation on superhorizon scales is well approximated in terms of that 
in the long wavelength limit, and that the latter can be constructed from the evolution 
of corresponding exactly homogeneous universe. Using these facts, we investigate the 
evolution of the cosmological perturbation on superhorizon scales in the universe dominated 
by oscillating multiple scalar fields which are generally interacting with each other, and the 
\Q • ratio of whose masses is incommensurable. Since the scalar fields oscillate rapidly around 
the local minimum of the potential, we use the action angle variables. We found that this 
problem can be formulated as the canonical perturbation theory in which the perturbed 
"q \ part appearing as the result of the expansion of the universe and the interaction of the 
scalar fields is bounded by the negative power of time. We show that by constructing the 
canonical transformations properly, the transformed hamiltonian becomes simple enough 
to be solved. As the result of the investigation using the long wavelength limit and the 
canonical perturbation theory, under the sufficiently general conditions, we prove that for 
the adibatic growing mode the Bardeen parameter stays constant and that for all the 
other modes the Bardeen parameter decays. From the viewpoint of the ergodic theory, it 
is discussed that as for the Bardeen parameter, the singularities appear probabilistically. 
This analysis serves the understanding of the evolution of the cosmological perturbations 
on superhorizon scales during reheating. 
PACS number(s):98.80.Cq 
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1 Introduction 



The cosmological structures such as the galaxies and the clusters of galaxy, etc, are thought 
to be produced from the seed pertubations through the growth induced by the gravitational 
instability. In the inflationary paradigm, the quantum fluctuation of the scalar field driv- 
ing the inflation called the inflaton plays the role of the seed perturbation. In the slow- 
rolling phase, the seed perturbations produced in the Hubble horizon are streched into the 
superhorizon scales by exponential cosmic expansion, and stay out of horizon until the per- 
turbation comes into horizon in the Friedman regime. We can calculate the amplitude of 
the seed pertubations when the potential of the inflaton is given, and know the evolutionary 
behavior of the cosmological perturbation in the horizon in the Friedmann universe when 
the composition of the cosmic matter is given. Therefore in order to obtain the knowledge 
of the inflationary universe from the statistical informations of the galaxy distribution, 
we need to know how the amplitudes of pertubations on superhorizon scales grow from 
the first horizon crossing in the slow rolling phase until the second horizon crossing in the 
Friedmann regime. 

The evolution of the cosmological perturbation on the superhorizon scales are thought 
to be simple. Because the Bardeen parameter are conserved in a good accuracy, when the 
equation of state is regular, and the energy of the universe is dominated by a single compo- 
nent. [3 ] But between the inflationary regime and the Friedmann regime, reheating occurs. 
During reheating, since the plural scalar fields oscillate rapidly around the local minimum 
of the potential, and since radiation is generated by the energy transfer from the scalar field 
oscillation, the equation of state becomes singular periodically [7 ] , and the entropy modes 
can appear [8 ]. Therefore the approximate conservation of the Bardeen parameter is not 
self-evident, and we need to solve the evolution of the cosmological perturbations without 
using the conservation of the Bardeen parameter. In general inflationary scenarios, there 
exisists plural scalar fields. [12 ] So in this paper, we consider the universe dominated by 
multiple scalar fields. 

Recently the strong method of solving the evolution of cosmological perturbations on 
superhorizon scales is developed. The evolutionary behavior of the cosmological perturba- 
tions on superhorizon scales is described by that in the long wavelength limit in a good 
accuracy. [7 ] Further the long wavelength limit of the solution of the evolution equation 
of the cosmological perturbations is generated from the perturbation of solution of the 
evolution equation of the background homogeneous universe. [4 ], [5 ], [6 ]. By using the 
fact, we investigate the evolution of the cosmological perturbation on superhorizon scales 
in the universe dominated by the oscillating multiple scalar fields interacting with each 
other. 

In order to use the above method, we are required to solve the exactly homogeneous 
system accurately. During reheating, the scalar fields oscillate rapidly around the local min- 
imum of the potential. As a method of treating dynamical system with rapidly oscillating 
degrees of freedom, we use the action angle variables, the hamiltonian formulation and the 
canonical perturbation theory. [9 ] In the canonical perturbation theory, the hamiltonian 
is decomposed into an unperturbed part whose solution is already known and an pertur- 
bation part, and by constructing the change of dynamical variables by which pertubation 
part is decreased, the system is made simple enough to be solved. The dynamical degrees 
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of freedom are decomposed into the slowly varying amplitudes of oscillation called action 
variables which are involutive first integrals of the unperturbed part and the fast increas- 
ing phase of oscillation called angle variables. The angle variables are eliminated from the 
hamiltonian by the canonical transformation. By using the method of the long wavelength 
limit and the canonical perturbation theory, we investigate the evolutionay behavior of the 
Bardeen parameter. While we adopt the action angle variables as the dynamical variables 
in the viewpoint of treating the reheating, from the standpoint decomposing the adiabatic 
and the entropy perturbation, the different variables were proposed. [11 ] Also in other 
approach the multiple scalar fields system was investigated. [10 ] 

This paper is organized as follows. In chapter 2, using the result of the paper [5 
], we explain the way to construct a solution of the perturbation equation in the long 
wavelength limit from a perturbation of a solution of the exactly homogeneous universe. 
In chapter 3, we write the system of the multiple scalar fields in the homogeneous universe 
in the hamiltonian form and fix the model which we treat concretely. In chapter 4, we 
apply the canonical perturbation theory to our system. The scalar fields oscillation is 
decomposed into the action angle variables, the former of which and the latter of which 
represent the slowly varying amplitudes and fast increasing phases of the scalar fields 
oscillation, respectively. By using the canonical transformations reccurently, we decrease 
the power of time variable which the interaction parts contain and eliminate fast oscillation 
from the evolution equation and extract the slow secular motion. We give the evaluation 
of the generating functions of the canonical transformations by which the transformed 
system becomes simple enough to be solved. In chapter 5, we investigate the evolution 
of the multiple scalar fields in the homogeneous universe. We prove the fact that the 
action variables are perpetually stable. In chapter 6, we investigate the evolution of the 
cosmological perturbations in the long wavelength limit paying attention to that of the 
Bardeen parameter. The upper bound of evolution of the perturbation variables and the 
uniqueness of the growing mode are derived. As a result, it is concluded that the Bardeen 
parameter stays constant in a good accuracy when the masses of the scalar fields are 
incommensurable. In this case, the Bardeen parameter becomes singular probabilistically. 
So in terms of the ergodic theory, we discuss the distribution of the singularity of the 
Bardeen parameter. The chapter 7 is devoted to the summary and the conclusion. 

Throughout the paper the natural units c — h — 1 are adopted and 8nG is denoted as 
k 2 . Further the notations for the perturbation variables adopted in the article [1 ] are used 
and their definitions are sometimes omitted except for those newly defined in this paper 

§2 Review; Evolution of cosmological perturbation in 
the longwavelength limit 

The dynamical behavior of cosmological perturbations on superhorizon scales is well ap- 
proximated by that of cosmological perturbations in the long wavelength limit [7 ], and 
the evolution of cosmological perturbations in the k — > limit can be constructed from 
the perturbation of the exactly homogeneous solution [5 ]. In this section, we present 
the results of the paper [5 ] in case of the universe dominated by multi-component scalar 
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fields. Next we explain two way to construct the perturbation of the exactly homogeneous 
solution. 

We only consider perturbations on a spatially Hat(K = 0) Robertson- Walker universe 
throughout the paper. Hence the background metric is given by 

ds 2 = -dt 2 + a{t) 2 dx 2 . (2.1) 

We consider the universe dominated by multi-component scalar fields, and the energy- 
momentum tensor is given by 

T£ = V> • V v <t> - (V> • V A + 2U) . (2.2) 

In the long wavelength limit, the evolution of the gauge invariant variable representing 
the fluctuation of the scalar field 4>i i n the flat time slice are described as 



HJ 2V \H' 
W(X 1 ,X 2 ) := X 1 -X 2 -X 1 -X 2 , (2.3) 

where Xi is a combination of exactly homogeneous perturbation variables described as 

Xi = 0<j>i - — — . (2.4) 
si a 

In the k — limit, it is noticed that the equality 

a 3 H 2 (6 \ 

~u~ w {w x ) =const (2 ' 5) 

holds. 

Corresponding to the way how the homogeneous solution is written, we can construct 
the exactly homogeneous perturbation in the different two ways. One is the way how it 
is given as derivative of the homogeneous solution with respect to the solution constants. 
The homogeneous solution is described as 

& = &(t,C), a = a(t,C), (2.6) 

or 

& = &(a,C), (2.7) 

where C is the set of all the solution constants characterizing homogeneous solution com- 
pletely. Xi is expressed as 

_ mt 1 C) ^l da(t,C) 
Xi dC Ha dC 1 } 



in the above case, and is expressed as 



A, = (2.9) 
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in the below case. Noticing the functional relation a = a(t,C), it is seen that the above 
two expressions are exactly the same. In this paper, the second expression is used when the 
evolution of the cosmological perturbations in the long wavelength limit is investigated. 

On the other hand, the solution is often expressed as union of the manifold on which 
the first integral is constant 

F„(0,^, a) = constant (n = 1, • • •, 2N) (2.10) 

where is canonical momentum conjugate to <fi as defined later. In this case, the exactly 
homogeneous perturbation is produced as 

Xi = {F n (0,p ,a),0i} = — * . (2.11) 

§3 Evolution equations of corresponding exactly ho- 
mogeneous universe 

In this section, we formulate the evolution of the exactly homogeneous universe dominated 
by multiple scalar fields as a hamiltonian form. Since our system are derived from the 
general covariant theory, it becomes a constrained system with gauge degree of freedom. 
By choosing the function of the scale factor as the time variable, we fix the gauge. It is 
the important first step since the evolution of the cosmological perturbations in the long 
wavelength limit is generated from the exactly homogeneous perturbations. 

We consider the situation where the multiple scalar fields interact with each other 

S = ( dt( Pa a + Y.P^i ~ N °h)- (3-1) 
J i 

Here N is the lapse function and Ch is the Hamiltonian constraint written as 

c *=^(-£tkrt +1 2m?X4.+m), (3-2) 

where f2 is the comoving volume of the constant time hypersurface. The potential U (0) is 
decomposed as 

^) = iE^ + M# (3-3) 

z i 

where Ui nt is assumed to be a sum of monomials of <fi of degree not smaller than 3 as in 
most physical context. We fix the gauge by 

= ip = a -t, (3.4) 

and its consistency condition with time evolution determines N as 

dib k 2 N . . 

= i = ~ aa^ ~ < 3 ' 5 > 



4 



For notational simplicity we write a in stead of t henceforth. By using Dirac bracket 
formalism, the Hamiltonian constraint = Ch is changed into a strong equation, and p a 
can be removed away. We obtain a reduced action with no constraint 

S = J2p^d(f)i - h a da, (3.6) 
J i 

h ° = ^M^XX + W)) 1 ' 2 (3-7) 
We make all physical quantities non-dimensional by dividing them with typical scales as 
C = (-) 3/2 , (3.8) 
$ = A , (3.9) 

00 

P* = 3 ~a~7v (3-10) 
/ii = — , (3.11) 



m 

e = ^^o, (3.12) 
S 



Qalm <f)Q ' 
J i 

-\ E ^ + ^c 2 e ^ 2 ^ 2 + e^uum) 1/2 dc (3.13) 

From now on, we use t instead of ( without possibility of confusion. In this paper, 
the evolution of the cosmological perturbations in the long wavelength limit is generated 
from the derivative of the exactly homogeneous solution with respect to solution constant 
with the scale factor fixed. Since we adopt the function of the scale factor as the time 
variable, in order to generate the evolution equations followed by the exactly homogeneous 
perturbations, it is sufficient to take variations with respect to the canonical variables with 
the time variable fixed in the unperturbed canonical equations of motion. 



§4 Integration Method 

This section is the main part of the paper. Since the scalar fields oscillate around the local 
minimum of the potential, we use the action angle variables, the former of which and the 
latter of which represent the slowly varying amplitudes and the fast increasing phases of 
scalar fields oscillation, respectively. We found that our system can be formulated as the 
canonical perturbation theory the perturbed part of which depending on the angle variables 
and the negative power of the time appears as a result of the expansion of the universe 
and the interaction of the scalar fields. Reducing the evolution equation into a simpler 
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form whose perturbed part is sufficiently small by constructing appropriate coordinate 
transformations is known as the method of variation of constants mathematically. Since 
our system is described in the hamiltonian form, we use the canonical transformations. By 
constructing the canonical transformation properly, we show that the perturbed parts of 
the hamiltonian can be decreased so that the transformed hamiltonian can be solved easily. 
The starting hamiltonian is given by 

« - + + e ^f^ 12 - (41) 

By performing canonical transformation ($;, P^f) — > rf) induced by generating function 

as 

W{*, rf) = ]T y/2y/mtri% + £ V**4*? ~ £ i\f\ (4.2) 

i i "" i 

we introduce complex canonical variables as 

e = ^(v^m--±=Up <p j, (4.3) 

rf = -U-^^ + ~P^). (4.4) 
Then the transfomed hamiltonian is given by 
H' = H + dW 



dt 

= ;OEMY + ^i^t) 1/2 + ^(E^ 2 + Ee 2 )- (4-5) 

Since the interaction term U- mt is written as a sum of monomials of $i; 

••^^Wfl' (4 ' 6) 



of degree not smaller than 3, we assume that U[ llt is a sum of monomials in terms of 
C/t = (C/tirf/t) °f power not smaller than 3 with t- independent coefficients; 



t 2 



,u mt = \e £ c kl dn v 7 y+xt 2 £ c k dn v 7 y, (4.7) 



m O0O |fc|+ji|=3 t t \k\ + \l\>4: t t 

where A is a constant characterizing the strength of interaction. We introduce the action 
angle variables by 

f = vOV* (4.8) 

rf = -is/Ti e - i9i . (4.9) 

Since k ^ I under the condition \k\ + \l\ = 3, the first term of the right hand side K 3 ; 

K 3 = jk 3 (J,6), (4.10) 



where k% is 27r-periodic in 9 — (0 1 , • • 9 N ) and is of power 3/2 in terms of J, satisfies 

; /•27T /*27r 

...y o d N oK*(j,o,t) = Q. (4.ii) 



(27T) 

The second term if 4 is written as 

^4 = Af;ifc 4p+2 (J,0), (4.12) 

p=2 J 

where fc 4 p+ 2 is 27r-periodic in 6 — (9 1 , • ■ 9 N ) and is of power (p + 2)/2 in terms of J. The 
transformed hamiltonian H' is rewritten as 

H = -(^■J + K 3 + K i ) 1/2 -J2 — ^0 l 
6 i t 

= 2 -^-J) l ' 2 + H. (4.13) 



The interaction part H is decomposed as 



H = -(//• J + K 3 + K A y/ 2 --(fi- J) 1 / 2 -^ — sin 20* 
e e i t 

= + (4.14) 



where 



Bi = 7(^V A ' 3 -?T si ' 129 " < 4 ' 16 > 



Now we assume that A/e is at most of order unity. By transferring ^-dependent part of A[ 
to B[, we obtain as a starting hamiltonian 

H=- e (fi- Ji) 1/2 + A 1 (J 1 , t) + Sx( Jx, 0!, t), (4.17) 
satisfying (i) // • J\ is bounded as 

lAi-Jil ><ri, (4.18) 
for some positive U\ on the domain (Ji, 1/t) G -D'i(pi) where 

^(P) = iU \); \J ~ Jno\ < P, \]\ < 1 + p, |Im^| < p}. (4.19) 

where J is the initial value of J; J = — !)■ 

(ii) t 2 Ai is analytic in terms of J±, 1/t and is real for real J±, 1/t and the inequality 

|t 2 Ai| < M} 15 (4.20) 
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holds for some positive on the domain Di(pi), where 

D n (p) = {(J, 6, j-); \J - J n0 \ < P, \Im9\ < p, |1| < 1 + p, |Imi| < p}. (4.21) 

(hi) tB\ is analytic in terms of J±, 9i, 1/t and is real for real J±, 9\, 1/t and is 
27r-periodic in 9\ — (9{, ■ • •, 9±) and the inequality 



tB x \ < M 2 (1) (4.22) 



holds for some positive and satisfy 



/•27T /"27T 

/ •••/ d N 9 1 tB 1 (J 1 ,9 1 ,t)=0, (4.23) 



(27T) 



N 



on the domain Di(pi). 

It is noticed that the perturbed parts A±, B\ appearing as the result of the expansion 
of the universe and the interactions of the scalar fields depend on the time as negative 
power. Since we are interested in the evolution for sufficiently large t, this property plays 
an essential role in the following analysis. 

In reheating phase, scalar fields oscillate rapidly. Since rapid oscillation is diffcult to 
treat, we construct change of coordinates in order to eliminate fast phases called the angle 
variables from the transformed equations of perturbed motion. The canonical perturbation 
theory is effective for this purpose, and can separate the slow secular motions from the fast 
oscillations by performing the canonical transformations recurrently, when the frequencies 
of unperturbed motion are incommensurable. With noticing that we can eliminate the 
dependence on the angle variables from the hamiltonian, it is important to understand 
that we can increase the power of the inverse of the time 1/t which is contained by the 
perturbed part depending on the angle variables B since we are interested in the evolution 
for sufficiently large t. 

Propositionl (Theorem for the canonical transformations) 

Suppose that m > 1. Suppose that the hamiltonian 

H {m) = - e (p- J m ) 1/2 + A m (J m , t) + B m ( J m , 9 m , t) (4.24) 

satisfy the following conditions: 

(i) For some positive C and d, the inequality 

\k-p\>CjL (4.25) 

holds for every integral- valued vector k = (ki, ■ ■ /cat). 

(ii) p • J m is bounded as 

\p ■ J m \ > cr m , (4.26) 

for some positive a m on the domain (J m , 1/t) e D' m (p m ). 

(hi) t 2 A m (J m ,t) is analytic and real for real J m , 1/t and the inequality 

\t 2 A m \ < M[ m) (4.27) 
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holds for some positive M 1 m ' ) on the domain (J m ,l/t) G D' m (p m ). 

(iv) t m B m (J m , 9 m , t) is periodic in m = (6^, ■ ■ ■, 0ffl). It is analytic and real for real 
J m , 9 m , 1/t and the inequality 

\t m B m \ < e™-^^ (4.28) 
holds for some positive on the domain (J m ,0 m , 1/t) G D m (p m ) and satisfy 

(2V^ L '"I d 6 rnt m B m (J m , 9 m , t) = 0. (4.29) 
For an arbitrary positive S, there exists eo such that, for an arbitrary e satisfying 

< e < e , (4.30) 

there exists the canonical transformation induced by the generating function S m ( J m +i, m , t) 
satisfying the following condition; 

(v) t m S m ( J m +i, m , t) is periodic in 9 m = (O^, ■ ■ ■, 0ffl). It is analytic and real for real 
Jm+ij @m, 1/t and the inequality 

\t m S m \ < e m L { ™ ] (4.31) 

holds for some positive on the domain (J m +i,0 m , 1/t) G D m+ i(p m+ i), where p m+ \ = 
Pm ~ 5. 

The transformed system defined by 

dS m 

Jm = Jm+1 + r, a , (4.32) 
0U m 



dS„ 

dJm+l 



l m+ i = m + —^, (4.33) 



H {m+1) = H (m) + 

at 

2 

— -(fj, • Jm+i) 1 / 2 + A m+1 (J m+ i,t) + B m+1 (J m+1 , 9 m+ i,t) (4.35) 

has the following properties; 

(vi) /x • J m +i is bounded as 

I At • J m +i\ > cr m+ i, (4.36) 

for some positive a m +i on the domain (J m +i, 1/t) G D^ +1 (p m+ i) 

(vii) t 2 A m+ i(J m+ i,t) is analytic and real for real J m +i, 1/t and the inequality 

|* 2 An+i| < M x (m+1) (4.37) 

holds for some positive M[ m+1 ^ on the domain (J m +i, 1/t) G -D^ +1 (p m+ i). 

(viii) t m+1 J B m+ i(J m+1 ,6' m+1 ,t) is periodic in 9 m+l = (6[ m+1 \- • -,0^ +1) ). It is analytic 
and real for real J m+ i, 9 m +i, 1/t and the inequality 

\t m+1 B m+1 \ < e m M 2 (m+1) (4.38) 
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holds for some positive on the domain (J m+ i,0 m+ i, 1/t) G D m+ i(p m+ i) and satisfy 

j p27r p2ir 

Jq^n j ' ' ' Jq ^ ®rn+it m B m+ i(J m+ i, 9 m+ i, t) — 0. (4.39) 

Proof 

We consider the canonical transformation induced by the generating function given by 

S m (J m+1 ,9 m ,t) = J2 S k (J m+1 ,t)e lk - 9 ™. (4.40) 

fc^O 



B m is decomposed as 



B m (j m+1 ,e m ,t) = Y,h(J m+ i,ty k - 9 -. (4.4i) 

fc^O 



The transformed hamiltonian is 

H {m+l) = H (m) + ° 



2 

= -(/i ■ J m+ i)^ 2 + A m ( J m+ l, t) 

i 1 1 dS m (J m+ i,9 m ,t) 

+ i? 2 + i? 3 + dSm ^ Jm +^ dm ^ (4.42) 



where R±, • • • R 3 are defined as 

i2i = "0* • ^) V2 " V W /2 " - ? r • ^ m(J ^ 1>6>m>t) , (4.43) 

e e e (// • J m+ i) 1/2 d0 m 

i? 2 = A m (J m , t) - A m (J m+1 ,t), (4.44) 

-R3 = B m (J m , 6 m , t) — B m (J m+1 , 9 rn , t). (4.45) 

We determine the generating function S m so that the leading term of the oscillating part 
can be eliminated; 

1 1 dS m (J m+ i,8 m ,t) . , 

~e (,-J m+1 y^ W m = - B Mm +1 ,e m ,t). (4.46) 

By comparing the fourier components in the both side, we obtain 

S k (J m+1 ,t) = ie(fi ■ J m +i) 1/2 7 —r\ b k (Jm+i, t). (4.47) 

(fx ■ k) 

Since from (iv), on the domain (J m +i, 1/t) e D' m (p m ), the inequality 

|* ro & fc (^m+i,f)| < e m - 1 M 2 (m) exp(-|A;|p m _ 1 ) (4.48) 
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holds, we evaluate the upper bound of the fourier components St, 

\t m S k (J m+1 ,t)\ < e^-Jrno + N/ipjV^e^M^eM-lklPrn) 

M (m) H 1 

< € ™ {fjL . j m0 + Nw ji/2^ { _y_ ew{ _ mPm _ S )}. (4.49) 
where in proceeding from the first row to the second row we use the inequality 

for arbitrary positive \k\, d, 5. On the domain | J m+ i — J m o| < Pm; | im #m| < Pm — 25, P"^™ 
is bounded as 

|r^ m | < £|rS fc |exp(|A;||Ini0 m |) 

< e m (p • J m0 + N^p m f/ 2 -^{-) d - £ exp(-|fc|5) 

° e k^O 

< e m 4 m) , (4.51) 
where is defined as 

4"° = (p • J m o + N ^Y /2 ^(-/y4) N (4-52) 
In passing from the second row to the third row, we use the inequality 

£exp(-|fc|*) < [l + 2Eexp(-A;5)] iV 

fc^O fc>0 

< (~ 5 ) N , (4-53) 

where we assume that 5 < 1. On the domain | J m +i — J m o\ < Pm, |Im# m | < p m — 35, 
|l/i| < 1 + p m , t m dS m /d9 rn is bounded as 

um dS m (J m+ i, 6 m , t) ^ ^ e m jj ^ ^ ^ 



So the inequality 



it j I _ | dS m (J m+ i o, 6> m o, t — 1) m Z^ 

KmO — 4+1 o — I 777. 1 S e — (4.00J 

dfl m o o 



holds. Of course, the above two inequalities hold on the domain | J m+ i — J m o | < p m — 35, 
\lm9 m \ < p m - 35, \l/t\ < 1 + p m . Then 

\Jrn+l — Jm+l o| < \Jm+l ~ Jmo\ + \JmO ~ Jm+l o| 

Am) 

< Pm -35 + e m ^-. (4.56) 
5 
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When we assume that e satisfies 

j(m) 

e mi < tfj ( 457 ) 

5 

the inequality 

| J m+ i - J m+ i | < p m - 25 (4.58) 

holds. The inequality (4.51) holds on the domain (J m +i,6 rn , 1/t) e -D(p m — 25). 

Now for an arbitrarily fixed J m+ i on the domain \J m +i — J m +i o| < Pm — 5, we consider 
the mapping T/ m+1 defined by 

Tj m+1 ;9 m — > m+ i = # m + -^-z (J m+ i,0 m ,f) (4.59) 



For an arbitrary 6> m on the domain |Im# m | < p m — 25, the inequality 



|$m ~ $m+l| < 1 777 {Jm+l)Qm-,~k) 



-Am) 

< e m (l + p m r^- (4.60) 

holds. Therefore the mapping Tj m+1 is injective (so diffeomorphic) on the domain |Im0 m | < 
p m — 25 — 4e m (l + p m ) m I '5 . Tj m+1 maps an arbitrary 9 rn on the domain |Im# m | < p rn — 
25-4e m (l + p m ) m 4 m) /5to 6 m+l on the domain |Im0 m+1 | < p m - 25 - 3e m (1 + p m ) m L { ™ ] / 5 . 
So on the domain \ J m +i - J m +i o| < Pm - 5, |Im0 m+ i| < p m - 25 - 3e m (l + p m ) m / '5 , 
the inverse mapping TJ^ is defined and diffeomorphic and |ImT7j" +i # m+ i| = |Im0 m | < 

p m - 25 - 4e m (l + PmTL^ 15. From now on, we assume that e satisfies 

j{m) 

3e m (l + p m ) m ^-<5. (4.61) 

Then the canonical transformation (J m +i,0 m +i) — > (J m ,0 m ) induced by the generating 
function S m (J m+ i,6 m ,t) is defined and is diffeomorphic on the domain \ J m+ \ — J m +i o| < 
p m - 5, \Im6 m \ < p m - 35 and llmTj^^+il = |rm0 m | < p m - 35 + 5/3 < p rn - 25. 

Next on the domain \J m +i — J m +i o| < Pm — 25, |Im# m | < p m — 35, we evaluate the 
residual part R defined as 

R = R 1+ R 2 + R 3+ dS m(Jm+l,9m,t) ^ 

at 



R\ can be rewritten as 

Ri 



f J m 1 1 1 

Jj m+ i ' dJ ~e { {p ■ J) 1/2 _ (^ . J m+1 )V2 > 

rJm 1 /■ 1 
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We evaluate 

t (m)2 1 

2m D I ^ ,2m-l ,,2 Ar2/i , „ \2m^l - 1 



|t 2m i?i| < e /j, N (1 + p„ 



A2 3/2 



T (m)2 1 

< e> 2 iV 2 (l + p m ) 2m 1 



S 2 at 2 ' 



Therefore we obtain 

J ( m ) 2 1 

m+ln I «< m,,2Ar2/-| , „ \3m-l ^1 1 



^ air 



As for i? 2 , we obtain 



|r +2 i? 2 | = r+ 2 f Jm dAm ^ mt) dj t 

^ ru2 9A m {J m , t) 

< sup{|t — \}\t 

< e m 



m | 

m+l u, Jm 

2 dA m ( y J m , t) ^^ m dS m { J m +i, 9 m , t) 
dJ m ' d6 m 

M (m) L {rn) 



Q5 2 

Therefore we obtain 



M (m) j (m) 

r +1 ^ 2 |<e m (l + Pm ) Ml Li 



65 2 



As for R 3 , we obtain 

dB m (J m ,6 m ,t) 



t 2m R 3 \ = \t 2m r """'^'"-""^ 4 

< GUp r I ^^^tt^ C ^r>T- ; t) i j |^ m ( 9'S'm('^m+li ^m; t) 



d.J m ' d6 m 



< e 



A Am) j (m) 
2m-l IVI 2 L>1 

65 2 

(m) j (m) 



ft/TV 11 ) T V 



68 2 

Therefore we obtain 



Also we obtain 



dt 

So R is evaluated as 



]\/r M t ( m ) 

| r+ l^|< £ m (l + pmr -1 ^2^1 



C 7" ( m ) 

|r+ l^! < e m [{m _ 1)L (m) + (1 + pm) ^ ; 



t m+1 i?| < e m M, 
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where 

7" (m) 

+ (m-l)4 m) + (l + p m )^-, (4.72) 

on the domain (J m+ i, m +i, G D m+ i(p m — 25). We decompose i? into the fourier 
components as 

R(J m+1 ,9 m+1 ,t) = J2 R k(Jrn + i,t)e ike ™+\ (4.73) 

k 

As for A m+ i, B m+ i defined as 

A m +i(J m +i,t) = A m (J m+ i,t) + Ro(J m+ i,t), (4.74) 

B m+1 (J m+1 ,6 m+1 ,t) = Y.MJm + ut)e ike ™ + \ (4.75) 



fe^O 



we obtain 



\t 2 A m+1 \<M[ m+1 \ (4.76) 

\t m+1 B rn+1 \ < e m M 2 (m+1) , (4.77) 

where 

M (m+1) = M (m) + (1 + pm )--2 emM; (4>7g) 

M 2 (m+1) = 2M (4.79) 

on the domain (J m +i, m +i, 1/i) G D m+ i(p m — 25). We put p m = p m _i — 25. We complete 
the proof. 
Proof End 

The following proposition is used in evaluating the difference between the original per- 
turbation variables and the transformed perturbation variables. Since as shown later, the 
action variable perturbations and the angle variable perturbations are bounded by t°, t re- 
spectively, the difference between the original perturbation variables and the transformed 
perturbation variables is bounded by some positive constant. 

Proposition2 

Suppose that the generating functions S m are analytic and the inequalities 

\t rn S m \ < e m 4 m) , (4.80) 
\t m (^—) ai (^r 2 S m \ < e m Lt ] («! + a 2 < 2) (4.81) 

holds for some positive L^, L^. 
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As for m > 2, the inequalities 



\5J[(t) - 5JUt)\ < jCi?\5J m (t)\ + € -Cg>\56 m (t)\, (4.82) 

\S0\(t) - 59^)1 < € -CiT ] \SJm(t)\ + -Cif\50 m (t)\ (4.83) 

holds for some positive C[™\ C[™\ where 

\SJ m (t)\= m^\Sr m (t)l (4.84) 

l<i<N 

\59 m (t)\= max : (4-85) 

l<i<A r 

Proof 

By taking the variation of 



we obtain 

■M 



<^ - M^-il < -iVLr j (|5J m+ i| + |<W m |), (4.87) 



I'd - <l < -NLt\\5J m+1 \ - |M„,|). i 1.88) 



As for \8Q l m \, we obtain 



WJ < K-^C-il + l^- 



m+1 1 



which yield 



< — iVLr j (|5J m+ i| + |^ m |) + |^ m+ i|, (4.89) 

\S6 m \ < ^NLt y \5J m+1 \ + (1 + e^^Lf ) ")|^ m+1 |. (4.90) 
By using the above inequality, the inequalities (4.87), (4.88) are rewritten as 

\8Jln - SJ L + i\ < ^NLt Y "(\6J m+ i\ + \59 m+1 \) : (4.91) 

m n+1 - 5911 < e ^NLt Y "(\5J m+1 \ + \56 m+1 \). (4.92) 

By using the above inequalities, by the induction, it is shown that for an arbitrary natural 
number m > 2, the inequalities 

\5J[(t) - 5.r m (t)\ < € jC[^\5J m (t)\ + e -C[^\59 m (t)\, (4.93) 

\Se\(t) - < \c^\5J m {t)\ + -C^\59 m {t)\ (4.94) 

holds for some positive C[™^ , , , ■ 



Proof End 
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§5 Evolution of the exactly homogeneous universe 

In this section, we investigate the evolutionary behavior of the exactly homogeneous uni- 
verse dominated by scalar fields. In terms of action angle variables, the evolution of the 
scalar fields is characterized by the perpetual stability of the action variables which we 
formulate below. 

Propositions (Perpetual stability of the action variables) 

Suppose that the hamiltonian 

H = ;(~E^ + ^T^*? + Am nt ($)) 1/2 (5.1) 

i i 

which satisfy the following conditions; 

(i) For some positive C and d, the inequality 

1*^1 >^ ( 5 -2) 

holds for every integral- valued vector k — (k±, ■ ■ k^)- 
In terms of action angle variables introduced by 

[~2 i—l 

$i = J — y Ji-cos9i, (5.3) 
y /ij t 

P^ = -Jfy i Jj i tsin9 i , (5.4) 



(ii) L r i n t( ( l > ) is a sum of monomials in terms of $ of power not smaller than 3. The 
inequality 

|At 2 f/ int ($)| <eM (5.5) 

holds for some positive M on the domain (J, 9, l/t) E D(p) 

(iii) The inequality 

I^E^ + ^ 2 E^I = I/^^I>- M) 

i i 

holds for some positive o on the domain (J, l/t) G D'(p). 

There exists some positive e such that for an arbitrary e satisfying 

< e < e , (5.7) 

the inequality 

\J-J \<eC (5.8) 

holds for some positive C. 
Proof 

We consider the sequence of canonical transformations defined in the Propositionl. 

H^(J 1 ,9 1 ,t)^H^(J 2 ,9 2 ,t). (5.9) 
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The transformed hamiltonian is given by 



H {2) = -y ■ J 2 ) 1/2 + A 2 ( J 2 , t) + B 2 (J 2 , 6 2 , t), (5.10) 



where 



The inequality 



(27T) 



! p2'K p2'K 

7f Jo "Jo ^3t 2 5 2 (J 2 ^ 2 ,t) = 0. (5.11) 



\t 2 A 2 \ < M {2 \ \t 2 B 2 \ < eM {2) (5.12) 
is satisfied on the domain D 2 (p 2 ). The equation of J 2 is 



dJj _ _dB 2 
It ~ ~~d¥ 2 

For an arbitrary positive 5, on the domain D(p 2 — 5) 

, dB 2 ] 

de. 



\m~4o\ < Jto\ 



2 



J £ M< 2) 



for an arbitrary t > 1. 

For an arbitrary positive e satisfying 



(5.13) 



M 2 (2) , . 

< e-l-. (5.14) 



M (2) 

v^e^- < p 2 - 5 (5.15) 



5 

the inequality 

\J 2 - J 20 \ ^VNe 1 -^-. (5.16) 
o 

holds. 

We complete the proof. 
Proof End 

§6 Evolution of Cosmological Perturbations in the Long 
Wavelength Limit 

We consider the evolution of the cosmological perturbations in the long wavelength limit 
since it approximates the evolution of the cosmological perturbations on superhorizon scales 
well [7 ]. In this section, by using the method of the long wavelength limit, and the canon- 
ical perturbation theory, we investigated the evolutionary behavior of the cosmological 
perturbations, in particular paying attention to that of the Bardeen parameter Z. We 



17 



show that for the adiabatic growing mode Z stays constant, and for all the other modes Z 
decays. Since the Bardeen parameter becomes singular probabilistically when the masses 
of the scalar fields are incommensurable, we discuss the distribution of the singularities of 
Z in the time axis in terms of the ergodic theory. 

On superhorizon scales, the Bardeen parameter defined by 



aH 

K-—a g , (6.1) 



is conserved in a good accuracy when the equation of state is regular and the energy of 
the universe is dominated by a single component. On the other hand, in the universe 
dominated by multi-component scalar fields, the Bardeen parameter is rewritten as 

Z = (6.2) 

(0) 2 

whose conservation is not self-evident. By using Eq.(2.5), we can see that the contribution 
from the second term in the right hand side of Eq.(2.3) is given by 



Z oc 



/«*!. (6.3) 



So we concentrate the contribution from the first term in the right hand side of Eq.(2.3) 
henceforth. 

Z = -H^M. (6.4) 
(0) 2 

The contribution to the Bardeen parameter from the exactly homogeneous perturbation, 
Eq.(6.4), is rewritten in terms of the action angle variables as 

Z = \-U ■ J + Kf' 2 1 



3t^ ' EiViJK 1 -cos 20*) 

]T(i<5J l sin 29 l - J\l - cos29 l )59 t ). (6.5) 

i 

In order to understand the evolutionary behavior of the perturbation variables such 

as Z, we have to investigate that of 59\ 5J l . Since the oscillating part is troublesome 

for investigating the differential equations of perturbations, we eliminate the fast phases 
by the canonical transformations. We consider the sequence of canonical transformations 
defined in the Propositionsl. 

j u 6l ,t)^ H&(J 2 , 9 2 , t) % H^ 3 \J 3 , 9 3 , t). (6.6) 
The transformed hamiltonian is given by 

# (3) = \ ■ hf 12 + A 3 (J 3 , t) + B 3 (J 3 , 9 3 , t), (6.7) 

where 

\t 2 A 3 \ < M{ 3) , \t s B 3 \ < e 2 M 2 (3) (6.8) 
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is satisfied on the domain D 3 {p 3 ). The variational equation of motion changed into the 
integral form is 



5Ji = 5X 



30 



59% 59: 



30 



t 



t 

+- I dt 



IfM f 
' t e Ji 



d 2 B 3 



d 2 B 3 



50 4- 
89 l 3 d9 3 3 d9\dJ 3 



dt5 



d 2 A 3 



t Ji dJldJ? 



6Js + \ldt[ 



d 2 B 3 

dJidO- 



-59 3 + 



3 2 B 3 



dJidJ; 



For an arbitrary positive 5, on the domain D 3 (p 3 — 5), we obtain 

\54(t)\ < 1*4,1 + |e 2 M 2 (3) {|^| m (t) + \\5JsUt)}, 



< m 

1 2 



+ ^^I^UW + ~M®\6J 3 \ m (t) 



e 2 a 3/2 ' 



+ -^e 2 Mi 3) {r-f\ m (t) + -\5J 3 \ m (t)}, 



'3 

,59 



1 



where 



They yield 



\A\ m (t) = sup \A(s)\ 

Ks<t 



\SJ 3 \ m (t) < (1 + eCOI^I + e 2 C 2 \59 30 \, 

ho I? 



|^U(t) < -C 3 |5J3o| + (l + eC 4 )|^ 3 



for some positive Ci, ■ ■ • C 4 . In particular, when Eq.(6.10) is written as 
59 i 3 = 59 i 30 + ^(t-l)V 3 (t) + W 3 i (t), 



Vs(t) = 



W, 



(t-1) 

^=1 



dtS 



dtJ^5J 3 + 



d.P 3 dJ 3 

from the inequalities (6.14), (6.15), we obtain 

1 



(/i-^3) 1 / 2 ' 



d 2 B 3 r „ d 2 B 3 



dJ l 3 d9 3 



59 3 + 



d.P 3 d.J 3 



\V 3 (t)-5- 



(/i- J30) 1 / 2 
m(t)\ <C 7 \5J 30 \+e 2 C 8 \59 



<e(C 5 |5J3o|+eC 6 |^3o| 



30| 



(6.9) 



(6.10) 
(6.11) 

(6.12) 
(6.13) 

(6.14) 
(6.15) 

(6.16) 
(6.17) 

(6.18) 

(6.19) 
(6.20) 



Next we have to know the difference between the original variables 59\, 5J\ and the 
transformed variables 59 3 , 5 J 3 , since we want to know the information of the original 
variables. As for the initial values between of 59±, 5J\ and of 59 m , 5J m , the inequalities 



m0 



<eCV\5Jio\ 



eC'^\59 w \ 



'10 



-59l n0 \<eCt ) \SJ: 



10 



+ eC /(m) ' 
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7 10 



(6.21) 
(6.22) 
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and the inequalities 

I^Jmol < (1 + eC^)\6J 10 \ + eC^\S9 10 \, (6.23) 

\59 m0 \ < eCif\5J 10 \ + (1 + eCi?)|^io|, (6.24) 

hold. Using (6.23), (6.24) to (6.14), (6.15) we obtain 

\SJ 3 \ m (t) < (1 + eC{ 4) )|5J 10 | + eCfl^iol, (6.25) 

I^U(i) < -Ci 4) \5J 10 \+Ci 4) \59 10 \. (6.26) 

By using (6.25), (6.26) to (4.82), (4.83), we obtain 

\5Jl(t) - 8Ji(t)\ < Ci|(JJio| + eC 2 |<f0i O |, (6.27) 

|^i(t) - S9i(t)\ < C 3 \5J 10 \ + eC 4 |^i |. (6.28) 

From (6.25), (6.27), we obtain 

\5Ji(t)\<C 5 \8J 10 \+eC 6 \89 10 \. (6.29) 

This yields the evaluation (6.37). 

On the other hand, from (6.19), (6.20), (6.23), (6.24), (6.28), we obtain the evaluation 
(6.37), (6.38), • • • (6.40). 

Proposition4 (Evaluation of the perturbation variables) 

Suppose that the hamiltonian 

H = -A~EPl + \ t2 H^M + Am nt ($)) 1/2 (6.30) 

i i 

which satisfy the following conditions; 

(i) For some positive C and d, the inequality 

\k^\>C^ d (6.31) 

holds for every integral- valued vector k — (k±, • • •, fcjv). 
In terms of action angle variables introduced by 

$i = J—yflXcoaOi, (6.32) 
y Hi t 

Pfc = -ij2fM'\fj i tsm0i, (6.33) 

(ii) C/i n t($) is a sum of monomials in terms of $ of power not smaller than 3. The 
inequality 

\\t 2 U hlt ^)\ < eM (6.34) 
holds for some positive M on the domain (J, 9, 1/t) e D(p). 
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(iii) The inequality 



I^E^ + ^ 2 E/^I = I/^I>- (6.35) 



holds for some positive a on the domain (J, 1/i) G D'(p) 
There exists some positive eo such that; 
for an arbitrary e satisfying 



< e < e , (6.36) 



the following inequalities holds. 
As for 5J l (t), the inequality 



\5J\t) - 5J l \ < d\5J \ + eC 2 \59 \ (6.37) 



holds for some positive C±, C 2 . 
When 59 % (i) is decomposed as 



60* (t) = 561 + !j(t - l)V(t) + W\t). (6.38) 



(i) V(t) is i-indendent and the inequality 

, 1 



\V(t) " ^r)l < <C 3 \5J \ + eC 4 \66 \) (6.39) 



holds for some positive C3 and C4. 

(ii) is i-dependent in general and the inequality 

\W\t)\<C 5 \5J \+eC e \59 \ (6.40) 

holds for some positive C5 and Cq . 

The decomposition of 59 % (t) into V(i) and W l (t) as in (6.38) is not unique. For the 
purpose explained below, it may be rather convenient to write in the next form. 

Proposition^ 

Under the same assumption presented in the above proposition, for t > 2, there exists 
a function V(t) such that 

8e\t) = 591 + ?±{t - l)V{t) + W'(t), (6.41) 

\V(t) -V(t = oo)\< y(Ci|<fJ | + eC 2 |^ |), (6.42) 

1^(01 < Csl^Jol + eC 4 |^ |, (6.43) 

for some positive Ci, • • •, C4. 
Proof 

We consider the transformed hamiltonian given by 

# (4) = hv J 4 ) 1/2 + ^4(^4, t) + B A { J 4 , 4 , *), ( 6 - 44 ) 
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where 

\t 2 A 4 \ < M{ 4) , \t 4 B 4 \ < e 3 M 2 (4) (6.45) 
is satisfied on the domain D 4 (p 4 ). As for J\{t) and 5Jl(t), we obtain inequalities 



roo 

\m- 4 j < j dt\ 



d6\ 



< e 3 C^, (6.46) 



roo f)^ R Pfl R 



e 2 



< ^(C 2 |5J | + eC 3 |^o|), (6.47) 
for some positive C±, C2, C3. As V(i) in the proposition, we adopt 

By using (6.46), (6.47), (6.48), for h,t 2 > 2, we obtain 

IVih) - V(t 2 )\ = \—^— t dt5- l —— + . ^"/^ 7T f 2 1 ' 



(/i-J^ 1 / 2 (*i - i)(f 2 - 1) A (/i-74) 1 / 2 ' 

+| (t 1 - 2 i)^ i -i ) r ^ (^3^^ • ^ 

"(*1-1)(*2-1) Z 2 ^ ( "2 ) (/.-J 4oo )3/ 2 ^ ' ^ ~' 

< Ir-rk'^l^ol+eCsl^ol), (6.49) 



which yields 



e 2 



- V(t = 00) I < -(C 4 \5J \ + eC 5 \56 \). (6.50) 

6 



We complete the proof. 
Proof End 

Using the above results to Eq.(6.5) we can see that for sufficiently large t, 

Z ^ --(fi- J + K) 1/2 V(t = 00). (6.51) 
o 

Therefore for the mode satisfying V(t = 00) 7^ which corresponds to the adiabatic 
growing mode, Z stays constant, and for other (2N — 1) modes satisfying V(t = 00) = 0, Z 
decays at most in proportion to 1/t. Then the uniqueness of the growing mode are derived. 
Because of the factor as 

£^J*(l-cos2^)' (6 ' 52) 
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all the terms in Eq.(6.5) except (6.51) become boundlessly large when the phase flow 9{t) 
approaches 9 l — or n (i — 1, N) boundlessly, although the numerator are not taken into 
account approaching zero. Does our system spend time near 9 l — or tt (i — 1, N) where 
the Bardeen parameter becomes large? The answer is as follows. Our track approaches 
the singular point boundlessly infinite times. The distribution of the singularities in the 
time axis is probabilistic and random. This problem can be regarded in the viewpoint of 
the ergodic theory. We formulate this as follows. 

Proposition5 

Under the condition 

(i)' If the equality /i ■ k — holds for some integral- valued vector k = (k 1 , ■ • •, k N ), k = 
holds. 

and the assumptions (ii), (iii) in the above proposition, 
there exists some positive e such that; 
for an arbitrary e satisfying 

< e < e , (6.53) 

the following statement holds. 

For all Riemannian integrable function f(9 1 ,---, 9 N ) which is 27r-periodic in 9 = (9 1 , ■ • 
•, 9 N ), the equality 

fenf rfW))=J ' (6 ' 54) 

holds where 

1= wp>C-C^ tm - (6 ' 55) 

We prove this proposition in the appendix. As a collorary of the proposition, we obtain 
Corollary 

Under the assumtions in the above proposition, as for a definition function <p(9) of the 
Jordan measurable set A C iJjlilP, 27r]; 

m = UlLvMW (6 - 56) 

the equality 

^ijh[ dtmt)) = i^r T{A) (6 ' 57) 

holds where t(A) is a measure of the set A. 
Proof 

When the set A is Jordan measurable, its definition function <p(9) is Riemannian inte- 
grable. Applying the above proposition to 4>{9), we obtain 

= ^r(A). (6.58) 

Proof End 
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From the above corollary, it can be seen that when the masses of the scalar fields are 
irrationally independent, each track distributes uniformly in the torus T N , more accurately 
saying, the period for which an arbitrary track stays in an arbitrary domain A is propor- 
tional to the measure of A. Let us adopt the domain which contains the singular point 
6 l — or 7r (i — 1, • • •, N) as A. Then we can evaluate how long the Bardeen parameter 
becomes singular. But we cannot say when the Bardeen parameter becomes singular. The 
occurrence of singularities is essentially probabilistic and random. We can find no rule in 
a sequence of times when Z becomes singular. 

§7 Summary and Discussion 

During the reheating, the multiple scalar fields oscillate arond the local minimum of the 
potential. In order to understand the evolutionary behavior of the cosmological perturba- 
tion on superhorizon scales during reheating, we investigate that in a universe dominated 
by multiple oscillating scalar fields. The dynamical behavior of the cosmological perturba- 
tions on superhorizon scales is well approximated in terms of that of the long wavelength 
limit, and the latter can be constructed from the evolution of corresponding exactly ho- 
mogeneous universe. These facts provide us with a strong tool to investigate the evolution 
of the cosmological perturbation on superhorizon scales including the conservation of the 
Bardeen parameter. In order to solve the evolutions of the cosmological perturbations in 
the long wavelength limit, we need know the evolution of the exactly homogeneous system. 
Since the scalar fields oscillate rapidly, we use the action angle variables. When this prob- 
lem is formulated in the hamiltonian form, as the result of the expansion of the universe 
and the interactions between the scalar fields the perturbed part depending on the angle 
variable appears. We found that it is bounded by the negative power of the time and that 
by choosing the canonical transformation properly, the perturbed part, in particular the 
power of the time, can be decreased so that the transformed hamiltonian can be solved 
easily. By using the long wavelength limit method, and the technique of the canonical 
perturbation theory, we prove the fact that if the ratio of frequncies of the scalar fields 
are incommensurable, for the adiabatic growing mode Bardeen parameter stays constant 
and for the other (2iV — 1) modes, the Bardeen parameter decays. In addition we show 
that the Bardeen parameter diverges probabilistically. From the viewpoint of the ergodic 
theory, we evaluate the period for which the Bardeen parameter diverges. 

In this paper, since the evolution of the cosmological perturbations are assumed to 
be generated from the perturbations of the exactly homogeneous solutions, we can only 
take into account the interactions between the homogenous modes, although in the real- 
istic reheating scenario the interactions between the non-homogeneous modes cannot be 
neglected and may be rather strong. When we think that the analysis of the latter are 
too complicated to be treated, the investigation of this paper are important, since we can 
investigate analytically, and does not specify the details of the interactions. Although the 
analysis are confined in the non-resonant case, the elimination of fast phases by the canon- 
ical transformation is thought to be still effective also in resonance case. The analysis in 
resonant case is in future publications. 

Our final goal is understanding of the evolutionary behaviors of the cosmological per- 
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turbations during reheating. We have shown the constancy of the Bardeen parameter in 
two extreme cases, that is, in the nearly integrable case corresponding to this paper, and 
the chaotic case where the discription by the two component fluid model [8 ] becomes good. 
So in order to achieve our purpose completely, general numerical calculations taking into 
account the interactions between the non-homogeneous modes, taking its stands between 
the above two extreme cases are necessary. 
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§8 Appendix; Mathematical Preliminary 

We give a little knowledge of the complex analysis and the fourier analysis, since it is 
used in the physical investigations of this paper. These propositions are fundamental when 
we evaluate the upper bound of the generating functions of the canonical transformations 
which decrease the size and the power of the time of the perturbed part of the hamiltonian. 

The following proposition is often used when we evaluate the derivative with respect to 
the canonical variables, for example the defference between the original variables and the 
transformed variables and the coefficients appearing in the perturbations of the canonical 
equations of motion. 

Proposition 

f(x) is analytic and the inequality 

\f(x)\<C (8.1) 

holds for some positive C on the domain D. Then for an arbitrary positive 5, in the domain 
D — 5 = {x E D; the 5 neighborhood of x is contained in D}, the inequalities 



, d 2 f , 2C 



e*<*T (8 ' 2) 



f(x) = ±-ldyl^-. (8.4) 



W^-P- (8 ' 3) 

hold. 

Proof 

Use the Cauchy equation 

ft rl 

2ni J y — x 

Proof End 

The following proposition guarantee the fact that the derivative with respect to the 
time decreases the power of the time. 
Proposition 

t m S(t) is analytic function of s — 1/t on the domain s E C(p) where 

C(p) = {s; \s\ < 1 + p, |Ims| < p} (8.5) 

and the inequality 

\t m S\ < C (8.6) 

holds for some positive C. Then for an arbitrary positive 5, on the domain C(p — 5) 
t m+1 dS/dt is analytic and the inequality 

K m+1 f \<C{m+ l -±P) (8.7) 

holds. 

Proof 
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By differentiating the equation 

—S(-)=fdz— -S(-), (8.8) 

s m V J z-sz m V v ' 

we obtain 

— *— S'(-) = -m—S(-) - / dz . 1 . 9 — g(-), (8.9) 

from which the proposition follows. 
Proof End 
Proposition 

F(Q) = J2k Fk exp ik ■ is periodic in — (0±, ■ ■ On) and analytic and the inequality 

|F(0)|<C (8.10) 
holds for some positive C on the domain 

|Im0|<p. (8.11) 

Then the fourier components F k satisfy 

\F k \ < C exp(-\k\p) (8.12) 

where \k\ is defined by 

\k\ = \h\ + ■ ■ ■ + \k N \. (8.13) 

Proof 

The fourier coefficients are given by 

Fk = j2^l "Jo d9 i--- d6 " e M-ik-0)F(0). (8.14) 

By changing the contour into 0i = Xi + ipi, (0 < x^ < 2tt, pi = — sgn(fcj)p) taking into 
account that contributions from the both sides cancel, F k is evaluated as 

< U e M-\h\p)C 

i 

= Cexp{-\k\p). (8.15) 

Proof End 

The following proposition guarantees that sufficiently small generating function defines 
the diffeomorphic canonical transformations. We omit the proof. 
Proposition 

Suppose that the mapping T;x~^yis defined on the domain D. When T is differ- 
entiable and detdy/dx ^ and \x — y\ < e on the domain D, T is diffeomorphic on the 
domain D — 4e. 
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§9 Appendix; Ergodic Property of Conditionally Pe- 
riodic Motion 



In chapter 6, we pointed out that the Bardeen parameter becomes singular at a certain point 
in the torus T N . In the viewpoint of the ergodic theory, we cannot say when the Bardeen 
parameter becomes singular, but can say how long the Bardeen parameter singular. In 
this appendix, we discuss the ergodic property as for the universe dominated by multiple 
oscillating scalar fields whose masses are incommensurable, i.e. irrationally independent. 
In our system, it can be shown that an arbitrary track is dense everywhere in the torus 
T N . It is referred to as a conditionary periodic motion. 

We consider the system proposed in Eq.(4.17). We assume that /i is a linearly indepen- 
dent vector on Z N ; 

fl .k = 0^k = 0. (9.1) 

As a first step, we prove the lemma below. 
Lemma 

lim — — f dte lk - m = 0, (9.2) 
t^oo T-lJi ' v ; 

Proof 

We consider the quantity as 

JdtF(J,9,t). (9.3) 

If G(J,9,t) satisfying 

i i dG „ , „. 

exists, the equality 

/ mJ , », t , = [GM «)]„ + / *[- ? §§ + ? §§ - f ] (9,, 

holds. When F = e ik ' e , G = e/i{fi ■ k) ■ (fj, ■ J) l / 2 e ik - 6 . Then there exists positive M, N such 
that 

| ^ £ dte ik e 1 < 6 T ]_ i [2M + N In T\ -> (T — > oo). (9.6) 

We complete the proof. 
Proo/ End 
Lemma 

For a triangular polynomial 

Pm(0) = E Ae*', (9.7) 



\<M 



the equality 



lim — — / dtP M (6) = Pm (9i 
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holds. 

Proof 

The both sides of Eq.(9.8) is linear. Therefore from the above lemma, Eq.(9.8) holds. 
Proof End 

Lemma; Weierstrass' theorem 

Suppose that /(x) is continuous on ]lili[ a i5 h]- For an arbitrary positive e, there exists 
a polynomial -P(x) such that 



N 

|/(x)-P(x)|<e xeJIki,]. (9.9) 



i=l 



Proof 

We consider the integral; 



Since 



J n = [\l-v 2 ) n dv, (9.10) 

JO 

J* = /\l - v 2 )"^ (0 < 5 < 1). (9.11) 
J s 

Jn> £0-- V ) ndv = ^T[> ( 9 - 12 ) 

J* < (l-5 2 r(l-5)<(l-5 2 r, (9.13) 

^<(n+l)(l-5 2 )", (9.14) 

we obtain 

lim = 0. (9.15) 

Without loss of generality, we can put < Oj < 6j < 1. We insert ctj, such that 
< oti < cii < bi < fa < 1. As a polynomial, we adopt 

Pn(x) = * jT /(u)[l - ( Wl - Xl ) 2 ] n ■ ■ • [1 - (u N - x N ) 2 ] n d Ul ■ ■ ■ du N , (9.16) 

where 

c = nKA]. (9.17) 

i=l 

P n (x) is a polynomial of 2n degree in terms of Xj, (i = 1, iV). We evaluate how accurately 
the polynomial P n (x) approximates the continuous function /(x). We introduce sets A, B 
defined by 

N 

A = l[[xi-6,Xi + 6] (0<6<1), (9.18) 
i=i 

N 

B = l[[a i ,b i ]- (9-19) 

i=i 
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First we evaluate 

i = jf /(u)[i - ( Ul - Xl ) 2 r ...[i-( UN - XN f} n d Ul ■■■du N 

-2 N (J n ff(x). (9.20) 
In terms of v = u — x, A is transferred into A' defined by 



N 



A' = nt-M- (9.21) 



1=1 

Therefore we obtain 



I = f /(x)[l - ( Ul - Xl ) 2 ] n ■ ■ ■ [1 - (u N - x N f] n d Ul ...du N 

J A' 

+ jf (/(v + x) - /(x))[l - ( Ul - a^) 2 ]" • • • [1 - (u N - x N f] n d Ul ■■■du N 

-2 N (J n ) N f(x). (9.22) 

Since /(x) is uniformly continuous on x G B, there exists positive 5 which depends only 
on e (0 < 5 < 1) such that 

|v| < ^5 (9.23) 

and 

|/(v + x)-/(x)| <e (xgB). (9.24) 

Therefore we obtain 

|I| < 2 JV {(J B ) JV - ( J„ - J:)^}M + e2*( Jn - J* n f, (9.25) 

where 

M = sup|/(x)|. (9.26) 
xec* 

Next we evaluate 

II = / /(u)[l - ( Ul - xx) 2 ] n • • • [1 - (u N - x N f] n d Ul ■ ■ ■ du N . (9.27) 

JC\A 

Since in terms of v = u — x, C, A are transferred into C", A' defined by 

N 

i=i 

N 

A' = n[-M. (9-28) 



and 



we obtain 



i=i 



C / C J D = nhl ) l] ! (9-29) 
i=i 

II| < M2 N {{J n ) N - (J n - J^}. (9.30) 
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From Eqs.(9.25), (9.30), we obtain 

|P B (x) - /(x)| < |I| + |II| < 2M{1 - (1 - J -f) N } + e(l - J -f) N . (9.31) 

From Eq.(9.15), the right hand side can be smaller than 2e when we adopt sufficiently large 
n. We complete the proof. 

Proof End 

Lemma 

Suppose that fiO 1 , • • -,9 N ) is continuous and 27r-periodic. For an arbitrary positive e, 
there exists a triangular polynomial P(9 1 , • • •, 9 N ); 

P(e\...,e N ) = ^+ £ a k e^\ (9.32) 

1 \k\<M 



such that 
Proof 

We define rf) by 



\f(9)-P(9)\<e. (9.33) 



where 



<P^,v) = Pi---PNf(9\--;9 N ), (9.34) 

& = A cos0*, (9.35) 
r] i = p i sin9 i . (9.36) 

When p? = Q + ^ 2 = 1, (i = 1, • • -,N), (/>(£, rf) coincide with f(9). Since 77) is 
continuous for £j, ^, from Weierstrass' theorem, there exists a polynomial for £, 77 P n (£, 77) 
which approximates uniformly 0(£, 77) in a rectangle containing Q + rjf = 1 (i — 1, • • •, JV). 
This completes the proof. 

Proof End 

Proposition 

For an arbitrary Riemannian integrable function f(9), the equality 

lim -J— F dtf{9(t))=J (9.37) 

T^oo i — 1 Ji 

holds. 

Proof 

For an arbitrary positive e, there exists triangular polynomials Pi, P2 such that 

Pi<f< P2, (9.38) 

P2-T\<e. (9.39) 

From Eq. (9.38), we obtain 

— ^ PMWt < I f(9(t))dt< ¥ ^jj i P2(9(t))dt. (9.40) 
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From the above inequality, by using the ergodic equality of the triangular polynomial, we 
obtain 

* S to^mf^f /(«(*))* 

< lim sup =1- f fW))dt< Pi (941) 

T ->oo T > T() I — 1 Jl 

which yields 

lim sup _J_ f T fm)d t - lim inf -J— T /(*(*))* < e (9.42) 
t — >ooy>y j — 1 Ji r —»oo t>t j — 1 Ji 

Since e is arbitrary, and P x < / < P 2 , we obtain 

lim -J— f T f{d{t))dt = J. (9.43) 

T^oo J — 1 Jl 

We complete the proof. 
Proof End 
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